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Abstract. This Letter reports on the photometric detection of transits of the Neptune-mass planet orbiting the 
nearby M-dwarf star GJ 436. It is by far the closest, smallest and least massive transiting planet detected so far. 
Its mass is slightly larger than Neptune's ai M — 22.6 ± 1.9 M®. The shape and depth of the transit lightcurves 
show that it is crossing the host star disc near its limb (impact parameter 0.84 ± 0.03) and that the planet size 
is comparable to that of Uranus and Neptune, R — 25200 ± 2200 km — 3.95 ± 0.35 i?®. Its main constituant is 
therefore very likely to be water ice. If the current planet structure models are correct, an outer layer of H/He 
constituting up to ten percent in mass is probably needed on top of the ice to account for the observed radius. 
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1. Introduction 

While over 200 extrasolar planets have been detected so 
far, the minority of them that transit their parent stars 
have the highest impact on our overall understanding of 
these objects (see review by Charbonneau et al. 2007). 
They are the only ones with accurate estimates of mass, 
radius, and, by inference, composition. Further precise 
monitoring of the brightest of these systems during pri- 
mary and secondary transits has even permitted the direct 
study of the planetary atmospheres (e.g. Charbonneau et 
al. 2002, Marley et al. 2007, Grillmair et al. 2007 ). Until 
now, this group was composed only of gaseous giant plan- 
plus the very massive HD 147506 b (Bakos et al. 
2007). 

The existence of smaller planets with masses of 5 - 
25 was recently uncovered by radial-velocity surveys 
(e.g. Butler et al. 2004; Udry et al. 2007), raising imme- 
diate questions about their constitution. Objects in this 
range of mass could be composed primarily of H and He 
gas, water ice, or refractory material (rock/iron). The de- 
tection of photometric transits of such an object would 
bring a preliminary answer to these questions and has 



thus been eagerly awaited until now. The photometric pre- 
cision needed to perform such a detection is beyond the 
capability of ground-based telescopes for solar-type stars 
and most of the planetary composition models. However, 
this is not the case for M dwarfs: their small radius 
makes possible the detection from the ground of transits 
of Neptune-sized, or even smaller, planets. Furthermore, 
existing radial- velocity surveys target relatively bright M- 
dwarfs, allowing small telescopes to be precise enough to 
carry out such a transit detection. In this context, we set 
up a photometric follow-up program of M-dwarfs known 
to harbor a close-in low mass planet. 

We report here the first result of our survey, the de- 
tection of transits of the Neptune-mass planet orbiting 
around the nearby M-dwarf star GJ 436 (Butler et al. 
2004, hereafter B04; Maness et al. 2007, hereafter M07) us- 
ing the 0.6m telescope at the Observatoire Frangois-Xavier 
Bagnoud (OFXB, Switzerland), and its confirmation with 
the 1.2m Euler telescope (La Silla) and the Im and 0.46m 
telescopes at the Wise Observatory (Israel). 

2. Observations and results 
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GJ 436 is a close {d = 10.2 pc) M2.5V star with V = 10.67. 
It has a low rotation velocity and does not exhibit par- 
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ticularly strong chromospheric activity nor photometric 
variability (B04), indicating an age greater than 3 Gyr. 
A periodic Doppler signal revealing the presence of a 
low mass planetary companion was reported by B04. The 
new Doppler measurements presented in M07 gave for the 
planet a minimum mass M sin i = 22.6 ± 1.9 M^, period, 
P = 2.64385 ± 0.00009 days, and indicated an eccentric 
orbit with e — 0.16 ±0.02. They also revealed a long term 
trend ('~ 1 m s^^ per year), indicating a possible distant 
companion. 

The star was photometrically monitored for transits of 
its close-in (a = 0.0285 AU) Neptune-mass planet in B04, 
and the authors concluded that complete transits across 
the star could be ruled out for gas giant compositions and 
should be considered as unlikely for solid compositions. 
We decided nonetheless to include this star in our target 
list and to observe it from OFXB, judging that the photo- 
metric light curve presented in B04 could not completely 
exclude a shallow and/or grazing transit. 

2.1. OFXB 0.6m telescope 

The OFXB is a small observatory located in St-Luc, 
Switzerland, mostly devoted to outreach activities. The 
telescope is a Newton 0.60m reflector providing a f/3.5 fo- 
cal ratio. The CCD camera is an Apogee AP47p, equipped 
with a Marconi 47-10 back-illuminated chip, providing a 
20'x20' field of view. This equipment has demonstrated 
its potential in exoplanet researches by its participation 
in the characterisation of the transiting planet WASP-2b 
(Cameron et al., 2007). GJ 436 was monitored during 8 
nights over 18 days between April 2nd and April 20th 2007 
for a total of 1108 useful exposures. We scheduled our ob- 
servations according to the transit windows expected from 
radial- velocity data, plus some short sequences at random 
phases to assess the photometric stability of the target. 
We observed in the y-band and defocused to reach 60s 
exposure time, a good trade-off between time sampling 
and scintillation mitigation (Gilliland & Brown, 1992). 

After a standard bias, dark and flatficld correction, all 
images were reduced with the IRAF/DAOPHOT aperture 
photometry software (Stetson, 1987), adapting the reduc- 
tion parameters to the FWHM of each image. Differential 
photometry was then performed using the flux of nearby 
stars for which a significant variability could be rejected. 
The rms of the OFXB photometry varies from 1.7 to 6 
mmag over the different nights. These changes are im- 
putable to differences in weather condition. No stellar vari- 
ability is seen in phase with the orbital period, a point 
already quoted in B04. Two clear transit-like events are 
present in April 2nd and 10th light curves, at the phase 
expected from the Doppler data. 

2.2. Wise Im and 0.46m telescopes 

To secure our tentative detection, wc observed GJ 436 at 
the Wise Observatory (Israel) on April 24th with two tele- 
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Fig. 1. OFXB (black) and Wise {red: Im, green: 46cm) 
photometry phase-folded using the cphemerids and period 
presented in Maness et al. (2007). 
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Fig. 2. Euler V-band transit photometry. The best-fit 
transit curve is superimposed in red. 



scopes - Im and 0.46m - simultaneously. We observed in 
the i?-band with the Im, and used no filter for the 0.46m. 
The reduction procedure was the same as above. Despite 
cloudy conditions, we managed to reach a precision high 
enough to clearly detect a transit egress at the expected 
phase with both instruments. 

Figure 1 presents the OFXB and Wise photometry 
phase-folded using the ephcmerides and period presented 
in M07. A clear transit-like event is visible at phase ~ 
0.007. 



2.3. Euler 1.2m telescope 

OFXB and Wise data together confirm the transiting sta- 
tus of GJ 436 b, but do not have high enough quality to 
firmly constrain the transit parameters. We gathered fur- 
ther observations from the Euler 1.2m telescope located 
at La Silla Observatory (Chile). Observations occured in 
photometric conditions on May 2nd during 5 hours, en- 
compassing the whole transit window at high airmass (1.8 
- 2.1). The same strategy as the one used at OFXB has 
been applied at Euler during this night (V^-band filter, 80s 
exposure time, defocus to 9" ) resulting in very accurate 
photometric time series, as can be seen in Fig. 2. The 
reduction procedure was the same than above. The rms 
outside the transit is ^1.2 mmag, while the expected de- 
viation taking into account scintillation and photon noise 
from the target and the reference stars is ~1 mmag. 
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3. Parameters of the planet GJ 436 b 

The lightcurves clearly indicate an almost grazing transit, 
since its duration is about two times shorter than for a 
central transit in front of a M2.5V star. The flux drop 
during the transit is only 0.6%. These two factors probably 
explain the non-detection by B04 (these authors state that 
their data exclude a central transit deeper than 0.4%). 

Since the Euler lightcurve is of a much superior quality 
than the others, and covers the whole transit, we use only 
these data for the determination of the parameters. We 
fitted a transit profile to the Euler data using the Mandel 
& Agol (2002) algorithm, the orbital elements in M07 and 
the quadratic limb darkening coefficients of Claret (2000) 
for Teff = 3500 K, log 5 = 4.5 and [Fe/H] = 0.0 in the V fil- 
ter. The mass of the parent star was adopted as M = 0.44 
± 0.04 Mq (sec discussion in M07). For main-sequence 
field M dwarfs of such low mass, the mass-radius relation 
is very tight. Observational constraints from M dwarfs in- 
tcrfcrometry (see Ribas ct al. 2006 and references therein) 
indicate M/Mq ~ R/Rq to within a few percent, and we 
use this mass-radius relation to set the primary radius. 
The models of Baraffe et al. (1998) for ages between 1 
and 10 Gyr would indicate radius values 0.02 Rq lower. 
On the other hand, the application of the radius calibra- 
tion from infrared luminosity and temperature by Kervella 
et al. (2004) indicates a higher value near R = 0.50 Rq. 
We thus adopt M = R = 0.44 ± 0.04 in solar units for 
the lightcurve fit. The remaining free parameters are the 
transit central epoch, the radius ratio and the orbital in- 
clination. 

The resiilts of the fit are given in Table 1 . The planet 
crosses the host star disc near its limb (impact parameter 
0.84 ± 0.03). The determination of the orbital inclination 
lifts the sini degeneracy on the planet's mass from the ra- 
dial velocity orbit, so that Mpi = 22.6 ± 1.9 Mq. Our best 
solution gives a radius of Rpi = 25200 ± 2200 km (3.95 
± 0.35 Rq) for the planet. The uncertainty is mainly due 
to that on the mass and radius of the primary. The for- 
mal uncertainties due to the photon noise of the lightcurve 
are very small. However, correlated systematics can cause 
much larger errors on the transit parameters (e.g. Pont, 
Zucker & Queloz 2006). Changing the reference stars in 
the Euler photometry, or using the lower-accuracy data 
from St-Luc and Wise, leads to changes of ~5% in the 
radius ratio. 

If the stellar radius is left as a free parameter in the 
lightcurve fit, the best-fit values are R = 0.46 Rq, Rpi = 
26500 km and i = 86°. This is an independent indication 
that the radius determination of the primary is basically 
correct. 

4. Discussion 

The measured radius of GJ 436 b is comparable to that 
of Neptune and Uranus. Figure 3 places it in the context 
of the mass-radius diagram for Solar System planets and 
transiting exoplanets. In this part of the mass-radius di- 
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Table 1. Parameters for the GJ 436 system, host star and 
transiting planet. *: from M07. 



agram, the position of a planet is a direct indication of 
its overall composition, while other factors such as tem- 
perature play only a minor role (see e.g. Fortney et al. 
2007). In the current paradigm, intermediate-mass plan- 
ets are composed of some or all of these four layers: an 
iron/nickel core, a silicate layer, an ice layer (H2O, CH4, 
NH3), and an H/He envelope. The mass and radius that 
we measure for GJ 436 b indicate that it is mainly com- 
posed of water ice. It is an "ice giant" planet like Uranus 
and Neptune rather than a small-mass gas giant or a very 
heavy "super-Earth" . It must have formed at a larger or- 
bital distance, beyond the "snow line" where the proto- 
planetary disc is cool enough for water to condensate, and 
subsequently migrated inwards to its present orbit. 

The temperature profile inside the planet is not ex- 
pected to modify this qualitative picture. The atmosphere 
of GJ 436 b must be hot: the equilibrium temperature is 
520 K to 620 K depending on the albedo, and a greenhouse 
effect may heat it to much larger temperatures. Tidal ef- 
fects from its eccentric orbit must also inject energy in its 
interior, but the iron, rock and water equations of state 
are not very sensitive to temperature at high pressure. 

We can ponder whether the planet has an H/He en- 
velope like the ice giants in the Solar System, or if its at- 
mosphere is composed mainly of water vapor. Our best-fit 
radius value places it slightly above the "pure ice" com- 
position mass-radius line of Fortney et al. (2007). A small 
H/He envelope may thus be required - even more if an 
iron/rock core is present as expected. At the upper end 
of the radius error bar, the H/He envelope would have to 
represent up to 10% in mass according to the models of 
Fortney et al. 2007 (see Figure 3). The lower end of the 
range is close to the mass-radius line for pure ice plan- 
ets. Water ice mixed with methane and ammonia is less 
dense than pure ice under high pressures, so the presence 
of a significant amount of these compounds within the ice 
could make the planet large enough despite a rock/iron 
core to account for the observed radius without invok- 
ing an H/He envelope. GJ 436 b could therefore be an 
"Ocean Planet" (Leger et al. 2004). Because of the high 
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Fig. 3. Planetary mass-radius diagram (adapted from 
Fortney et al. 2007) comparing the position of Solar 
System planets, transiting hot Jupiters (diamonds), and 
GJ 436 b. The lines indicate the position of the Fortney 
et al. models for different compositions: pure iron, pure 
silicate, pure water ice (with thermal profiles from Solar 
System planets), and models for irradiated planets at 0.1 
AU from a Solar- type star with a fraction of 10%, 50% 
and 100% of Hydrogen/Helium. The dotted lines show 
the models for a cold (a = 10 AU) and very hot (a = 0.02 
AU) pure H/He gas giant. 



surface temperature, this would imply a steam atmosphere 
above supercritical water rather than an Earth-like situ- 
ation. As methane and ammonia have very low conden- 
sation temperatures, this scenario would imply migration 
from a wide orbit. 

It could be expected that on such a close orbit, an 
H/He envelope would quickly evaporate. But although the 
planet is very close to its parent star, the small size and 
low temperature of the primary mean that such an en- 
velope could be retained over long timescales (see M07, 
Lecavelier 2007). A more precise radius determination can 
help determine whether the planet has a water or H/He 
envelope. 



The fact that the orbit of GJ 436 b is not circular indi- 
cates a high tidal quality factor Q for the planet, compat- 
ible with an ice giant rather than a predominantly rocky 
planet - although the eccentricity could be due to the in- 
fluence of an unseen planetary companion, as pointed out 
by M07. 

GJ 436 b is the first hot Neptune with a radius mea- 
surement, and turns out to be a Neptune-like ice giant, 
mostly composed of water ice, not a rock/iron "super- 
Earth" , nor a low-mass gas giant. Its detection illus- 
trates the potential of extensive high-precision photomet- 
ric follow-up of planets detected by radial velocity. It is the 
closest transiting planet known, and opens many oppor- 
tunities for further observations to characterize the planet 
itself. 
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